C. elegans as a model for inter-individual variation in metabolism.

Individuals can exhibit differences in metabolism that are caused by the interplay of genetic background, nutritional input, microbiota and other environmental factors1-4. It is difficult to connect differences in metabolism to genomic variation and derive underlying molecular mechanisms in humans, owing to differences in diet and lifestyle, among others. Here we use the nematode Caenorhabditis elegans as a model to study inter-individual variation in metabolism. By comparing three wild strains and the commonly used N2 laboratory strain, we find differences in the abundances of both known metabolites and those that have not to our knowledge been previously described. The latter metabolites include conjugates between 3-hydroxypropionate (3HP) and several amino acids (3HP-AAs), which are much higher in abundance in one of the wild strains. 3HP is an intermediate in the propionate shunt pathway, which is activated when flux through the canonical, vitamin-B12-dependent propionate breakdown pathway is perturbed5. We show that increased accumulation of 3HP-AAs is caused by genetic variation in HPHD-1, for which 3HP is a substrate. Our results suggest that the production of 3HP-AAs represents a 'shunt-within-a-shunt' pathway to accommodate a reduction-of-function allele in hphd-1. This study provides a step towards the development of metabolic network models that capture individual-specific differences of metabolism and more closely represent the diversity that is found in entire species.

Screening by deep sequencing reveals mediators of microRNA tailing in C. elegans.

microRNAs are frequently modified by addition of untemplated nucleotides to the 3' end, but the role of this tailing is often unclear. Here we characterize the prevalence and functional consequences of microRNA tailing in vivo, using Caenorhabditis elegans. MicroRNA tailing in C. elegans consists mostly of mono-uridylation of mature microRNA species, with rarer mono-adenylation which is likely added to microRNA precursors. Through a targeted RNAi screen, we discover that the TUT4/TUT7 gene family member CID-1/CDE-1/PUP-1 is required for uridylation, whereas the GLD2 gene family member F31C3.2-here named GLD-2-related 2 (GLDR-2)-is required for adenylation. Thus, the TUT4/TUT7 and GLD2 gene families have broadly conserved roles in miRNA modification. We specifically examine the role of tailing in microRNA turnover. We determine half-lives of microRNAs after acute inactivation of microRNA biogenesis, revealing that half-lives are generally long (median = 20.7 h), as observed in other systems. Although we observe that the proportion of tailed species increases over time after biogenesis, disrupting tailing does not alter microRNA decay. Thus, tailing is not a global regulator of decay in C. elegans. Nonetheless, by identifying the responsible enzymes, this study lays the groundwork to explore whether tailing plays more specialized context- or miRNA-specific regulatory roles.

Paracellular and Transcellular Leukocytes Diapedesis Are Divergent but Interconnected Evolutionary Events.

Infiltration of the endothelial layer of the blood-brain barrier by leukocytes plays a critical role in health and disease. When passing through the endothelial layer during the diapedesis process lymphocytes can either follow a paracellular route or a transcellular one. There is a debate whether these two processes constitute one mechanism, or they form two evolutionary distinct migration pathways. We used artificial intelligence, phylogenetic analysis, HH search, ancestor sequence reconstruction to investigate further this intriguing question. We found that the two systems share several ancient components, such as RhoA protein that plays a critical role in controlling actin movement in both mechanisms. However, some of the key components differ between these two transmigration processes. CAV1 genes emerged during Trichoplax adhaerens, and it was only reported in transcellular process. Paracellular process is dependent on PECAM1. PECAM1 emerged from FASL5 during Zebrafish divergence. Lastly, both systems employ late divergent genes such as ICAM1 and VECAM1. Taken together, our results suggest that these two systems constitute two different mechanical sensing mechanisms of immune cell infiltrations of the brain, yet these two systems are connected. We postulate that the mechanical properties of the cellular polarity is the main driving force determining the migration pathway. Our analysis indicates that both systems coevolved with immune cells, evolving to a higher level of complexity in association with the evolution of the immune system.

Using Convolutional Neural Networks to Measure the Physiological Age of Caenorhabditis elegans.

Caenorhabditis elegans (C. elegans) is a popular and excellent model for studies of aging due to its short lifespan. Methods for precisely measuring the physiological age of C. elegans are critically needed, especially for antiaging drug screening and genetic screening studies. The effects of various antiaging interventions on the rate of aging in the early stage of the aging process can be determined based on the quantification of physiological age. However, in general, the age of C. elegans is evaluated via human visual inspection of morphological changes based on personal experience and subjective judgment. For example, the rate of motor activity decay has been used to predict lifespan in early- to mid-stage aging. Using image processing, the physiological age of C. elegans can be measured and then classified into periods or classes from childhood to elderhood (e.g., 3 periods comprising days 0-2, 4-6 and 10-12) by using texture entropy (Shamir, L. et al., 2009). Our dataset consists of 913 microscopic images of C. elegans, with approximately 60 images per day from day 1 to day 14 of adulthood. We present quantitative methods to measure the physiological age of C. elegans with convolution neural networks (CNNs), which can measure age with a granularity of days rather than periods. The methods achieved a mean absolute error (MAE) of less than 1 day for the measured age of C. elegans. In our experiments, we found that after training and testing our dataset, 5 popular CNN models, 50-layer residual network (ResNet50), InceptionV3, InceptionResNetV2, 16-layer Visual Geometry Group network (VGG16) and MobileNet, measured the physiological age of C. elegans with an average testing MAE of 1.58 days. Furthermore, based on the results, we propose two models, one model for linear regression analysis and the other model for logistic regression, that combine a CNN model and a new attribute: curved_or_straight. The linear regression analysis model achieved a test MAE of 0.94 days; the logistic regression model achieved an accuracy of 84.78 percent with an error tolerance of 1 day.

Functionally non-redundant paralogs spe-47 and spe-50 encode FB-MO associated proteins and interact with him-8.

The activation of C. elegans spermatids to crawling spermatozoa is affected by a number of genes including spe-47. Here, we investigate a paralog to spe-47: spe-50, which has a highly conserved sequence and expression, but which is not functionally redundant to spe-47. Phylogenetic analysis indicates that the duplication event that produced the paralogs occurred prior to the radiation of the Caenorhabditis species included in the analysis, allowing a long period for the paralogs to diverge in function. Furthermore, we observed that knockout mutations in both genes, either alone or together, have little effect on sperm function. However, hermaphrodites harboring both knockout mutations combined with a third mutation in the him-8 gene are nearly self-sterile due to a sperm defect, even though they have numerous apparently normal sperm within their spermathecae. We suggest that the sperm in these triple mutants are defective in fusing with oocytes, and that the effect of the him-8 mutation is unclear but likely due to its direct or indirect effect on local chromatin structure and function.

Whole-organism behavioral profiling reveals a role for dopamine in state-dependent motor program coupling in C. elegans.

Animal behaviors are commonly organized into long-lasting states that coordinately impact the generation of diverse motor outputs such as feeding, locomotion, and grooming. However, the neural mechanisms that coordinate these distinct motor programs remain poorly understood. Here, we examine how the distinct motor programs of the nematode C. elegans are coupled together across behavioral states. We describe a new imaging platform that permits automated, simultaneous quantification of each of the main C. elegans motor programs over hours or days. Analysis of these whole-organism behavioral profiles shows that the motor programs coordinately change as animals switch behavioral states. Utilizing genetics, optogenetics, and calcium imaging, we identify a new role for dopamine in coupling locomotion and egg-laying together across states. These results provide new insights into how the diverse motor programs throughout an organism are coordinated and suggest that neuromodulators like dopamine can couple motor circuits together in a state-dependent manner.

Animals generate many different motor programs (such as moving, feeding and grooming) that they can alter in response to internal needs and environmental cues. These motor programs are controlled by dedicated brain circuits that act on specific muscle groups. However, little is known about how organisms coordinate these different motor programs to ensure that their resulting behavior is coherent and appropriate to the situation. This is difficult to investigate in large organisms with complex nervous systems, but with 302 brain cells that control 143 muscle cells, the small worm Caenorhabditis elegans provides a good system to examine this question. Here, Cermak, Yu, Clark et al. devised imaging methods to record each type of motor program in C. elegans worms over long time periods, while also dissecting the underlying neural mechanisms that coordinate these motor programs. This constitutes one of the first efforts to capture and quantify all the behavioral outputs of an entire organism at once. The experiments also showed that dopamine ­ a messenger molecule in the brain ­ links the neural circuits that control two motor programs: movement and egg-laying. A specific type of high-speed movement activates brain cells that release dopamine, which then transmits this information to the egg-laying circuit. This means that worms lay most of their eggs whilst traveling at high speed through a food source, so that their progeny can be distributed across a nutritive environment. This work opens up the possibility to study how behaviors are coordinated at the level of the whole organism ­ a departure from the traditional way of focusing on how specific neural circuits generate specific behaviors. Ultimately, it will also be interesting to look at the role of dopamine in behavior coordination in a wide range of animals.

Brn3/POU-IV-type POU homeobox genes-Paradigmatic regulators of neuronal identity across phylogeny.

One approach to understand the construction of complex systems is to investigate whether there are simple design principles that are commonly used in building such a system. In the context of nervous system development, one may ask whether the generation of its highly diverse sets of constituents, that is, distinct neuronal cell types, relies on genetic mechanisms that share specific common features. Specifically, are there common patterns in the function of regulatory genes across different neuron types and are those regulatory mechanisms not only used in different parts of one nervous system, but are they conserved across animal phylogeny? We address these questions here by focusing on one specific, highly conserved and well-studied regulatory factor, the POU homeodomain transcription factor UNC-86. Work over the last 30 years has revealed a common and paradigmatic theme of unc-86 function throughout most of the neuron types in which Caenorhabditis elegans unc-86 is expressed. Apart from its role in preventing lineage reiterations during development, UNC-86 operates in combination with distinct partner proteins to initiate and maintain terminal differentiation programs, by coregulating a vast array of functionally distinct identity determinants of specific neuron types. Mouse orthologs of unc-86, the Brn3 genes, have been shown to fulfill a similar function in initiating and maintaining neuronal identity in specific parts of the mouse brain and similar functions appear to be carried out by the sole Drosophila ortholog, Acj6. The terminal selector function of UNC-86 in many different neuron types provides a paradigm for neuronal identity regulation across phylogeny. This article is categorized under: Gene Expression and Transcriptional Hierarchies > Regulatory Mechanisms Invertebrate Organogenesis > Worms Nervous System Development > Vertebrates: Regional Development.

Deep sampling of Hawaiian Caenorhabditis elegans reveals high genetic diversity and admixture with global populations.

Hawaiian isolates of the nematode species Caenorhabditis elegans have long been known to harbor genetic diversity greater than the rest of the worldwide population, but this observation was supported by only a small number of wild strains. To better characterize the niche and genetic diversity of Hawaiian C. elegans and other Caenorhabditis species, we sampled different substrates and niches across the Hawaiian islands. We identified hundreds of new Caenorhabditis strains from known species and a new species, Caenorhabditis oiwi. Hawaiian C. elegans are found in cooler climates at high elevations but are not associated with any specific substrate, as compared to other Caenorhabditis species. Surprisingly, admixture analysis revealed evidence of shared ancestry between some Hawaiian and non-Hawaiian C. elegans strains. We suggest that the deep diversity we observed in Hawaii might represent patterns of ancestral genetic diversity in the C. elegans species before human influence.

Origination and evolution of orphan genes and de novo genes in the genome of Caenorhabditis elegans.

Orphan genes that lack detectable homologues in other lineages could contribute to a variety of biological functions. However, their origination and function mechanisms remain largely unknown. Herein, through a comprehensive and systematic computational pipeline, we identified 893 orphan genes in the lineage of C. elegans, of which only a low fraction (0.9%) were derived from transposon elements. Six new protein-coding genes that de novo originated from non-coding DNA sequences in the genome of C. elegans were also identified. The authenticity and functionality of these orphan genes and de novo genes are supported by three lines of evidences, consisting of transcriptional data, and in silico proteomic data, and the fixation status data in wild populations. Orphan genes and de novo genes exhibited simple gene structures, such as, short in protein length, of fewer exons, and are frequently X-linked. RNA-seq data analysis showed these orphan genes are enriched with expression in embryo development and gonad, and their potential function in early development was further supported by gene ontology enrichment analysis results. Meanwhile, de novo genes were found to be with significant expression in gonad, and functional enrichment analysis of the co-expression genes of these de novo genes suggested they may be functionally involved in signaling transduction pathway and metabolism process. Our results presented the first systematic evidence on the evolution of orphan genes and de novo origin of genes in nematodes and their impacts on the functional and phenotypic evolution, and thus could shed new light on our appreciation of the importance of these new genes.

The Emergence of the Tardigrade Hypsibius exemplaris as a Model System.

The success of scientists in revealing biological mechanisms has depended in large part on choosing tractable model systems. In 1997, molecular phylogenetics revealed that two of biology's most tractable models-Caenorhabditis elegans and Drosophila-are much more closely related to each other than had been thought previously. I began to explore whether any of the little-studied members of this branch of the tree of life might serve as a new model for comparative biology that could make use of the rich and ongoing sources of information flowing from C. elegans and Drosophila research. Tardigrades, also known as water bears, make up a phylum of microscopic animals. The tardigrade Hypsibius exemplaris (recently disambiguated from a closely related species, Hypsibius dujardini) can be maintained in laboratories and has a generation time of <2 wk at room temperature. Stocks of animals can be stored frozen and revived. The animals and their embryos are optically clear, and embryos are laid in groups, with each synchronous clutch of embryos laid in a clear molt. We have developed techniques for laboratory study of this system, including methods for microinjection of animals, immunolocalization, in situ hybridization, RNA interference, transcriptomics, and methods for identifying proteins that mediate tolerance to extreme environments. Here, I review the development of this animal as an emerging model system, as well as recent molecular studies aimed at understanding the evolution of developmental mechanisms that underpin the evolution of animal form and at understanding how biological materials can survive extreme environments.

Improved strategy for the curation and classification of kinases, with broad applicability to other eukaryotic protein groups.

Despite the substantial amount of genomic and transcriptomic data available for a wide range of eukaryotic organisms, most genomes are still in a draft state and can have inaccurate gene predictions. To gain a sound understanding of the biology of an organism, it is crucial that inferred protein sequences are accurately identified and annotated. However, this can be challenging to achieve, particularly for organisms such as parasitic worms (helminths), as most gene prediction approaches do not account for substantial phylogenetic divergence from model organisms, such as Caenorhabditis elegans and Drosophila melanogaster, whose genomes are well-curated. In this paper, we describe a bioinformatic strategy for the curation of gene families and subsequent annotation of encoded proteins. This strategy relies on pairwise gene curation between at least two closely related species using genomic and transcriptomic data sets, and is built on recent work on kinase complements of parasitic worms. Here, we discuss salient technical aspects of this strategy and its implications for the curation of protein families more generally.

Evolutionary plasticity of the NHL domain underlies distinct solutions to RNA recognition.

RNA-binding proteins regulate all aspects of RNA metabolism. Their association with RNA is mediated by RNA-binding domains, of which many remain uncharacterized. A recently reported example is the NHL domain, found in prominent regulators of cellular plasticity like the C. elegans LIN-41. Here we employ an integrative approach to dissect the RNA specificity of LIN-41. Using computational analysis, structural biology, and in vivo studies in worms and human cells, we find that a positively charged pocket, specific to the NHL domain of LIN-41 and its homologs (collectively LIN41), recognizes a stem-loop RNA element, whose shape determines the binding specificity. Surprisingly, the mechanism of RNA recognition by LIN41 is drastically different from that of its more distant relative, the fly Brat. Our phylogenetic analysis suggests that this reflects a rapid evolution of the domain, presenting an interesting example of a conserved protein fold that acquired completely different solutions to RNA recognition.

The multipotency-to-commitment transition in Caenorhabditis elegans-implications for reprogramming from cells to organs.

In animal embryos, cells transition from a multipotential state, with the capacity to adopt multiple fates, into an irreversible, committed state of differentiation. This multipotency-to-commitment transition (MCT) is evident from experiments in which cell fate is reprogrammed by transcription factors for cell type-specific differentiation, as has been observed extensively in Caenorhabditis elegans. Although factors that direct differentiation into each of the three germ layer types cannot generally reprogram cells after the MCT in this animal, transcription factors for endoderm development are able to do so in multiple differentiated cell types. In one case, these factors can redirect the development of an entire organ in the process of "transorganogenesis". Natural transdifferentiation also occurs in a small number of differentiated cells during normal C. elegans development. We review these reprogramming and transdifferentiation events, highlighting the cellular and developmental contexts in which they occur, and discuss common themes underlying direct cell lineage reprogramming. Although certain aspects may be unique to the model system, growing evidence suggests that some mechanisms are evolutionarily conserved and may shed light on cellular plasticity and disease in humans.

Rapid Gene Family Evolution of a Nematode Sperm Protein Despite Sequence Hyper-conservation.

Reproductive proteins are often observed to be the most rapidly evolving elements within eukaryotic genomes. The major sperm protein (MSP) is unique to the phylum Nematoda and is required for proper sperm locomotion and fertilization. Here, we annotate the MSP gene family and analyze their molecular evolution in 10 representative species across Nematoda. We show that MSPs are hyper-conserved across the phylum, having maintained an amino acid sequence identity of 83.5-97.7% for over 500 million years. This extremely slow rate of evolution makes MSPs some of the most highly conserved genes yet identified. However, at the gene family level, we show hyper-variability in both gene copy number and genomic position within species, suggesting rapid, lineage-specific gene family evolution. Additionally, we find evidence that extensive gene conversion contributes to the maintenance of sequence identity within chromosome-level clusters of MSP genes. Thus, while not conforming to the standard expectation for the evolution of reproductive proteins, our analysis of the molecular evolution of the MSP gene family is nonetheless consistent with the widely repeatable observation that reproductive proteins evolve rapidly, in this case in terms of the genomic properties of gene structure, copy number, and genomic organization. This unusual evolutionary pattern is likely generated by strong pleiotropic constraints acting on these genes at the sequence level, balanced against expansion at the level of the whole gene family.

Gain and Loss of Small RNA Classes-Characterization of Small RNAs in the Parasitic Nematode Family Strongyloididae.

The nematode family Strongyloididae is of particular interest because it contains important parasites of medical and veterinary relevance. In addition, species of this family can form parasitic and free-living generations and it also occupies an interesting phylogenetic position within the nematodes. Nematodes differ in several ways from other taxa with respect to their small noncoding RNAs. Recent comparative studies revealed that there is also considerable variability within the nematodes. However, no Strongyloididae species or close relative was included in these studies. We characterized the small RNAs of two developmental stages of three different Strongyloididae species and compared them with the well-studied free-living nematodes Caenorhabditis elegans and Pristionchus pacificus. Strongyloididae have conserved and taxon-specific microRNAs, many of which are differentially regulated between the two developmental stages. We identified a novel class of around 27-nucleotide-long RNAs starting with 5'G or A, of which a large fraction have the potential to target transposable elements. These RNAs most likely have triphosphates at their 5' ends and are therefore presumably synthesized by RNA-dependent RNA polymerases. In contrast to C. elegans but similarly to some other nematode taxa, Strongyloididae have no Piwi-interacting RNAs, nor do their genomes encode Argonaute proteins of the Piwi family. Finally, we attempted but failed to detect circulating parasite small RNAs in the blood of hosts.

Comparative Transcriptomics of Steinernema and Caenorhabditis Single Embryos Reveals Orthologous Gene Expression Convergence during Late Embryogenesis.

Cells express distinct sets of genes in a precise spatio-temporal manner during embryonic development. There is a wealth of information on the deterministic embryonic development of Caenorhabditis elegans, but much less is known about embryonic development in nematodes from other taxa, especially at the molecular level. We are interested in insect pathogenic nematodes from the genus Steinernema as models of parasitism and symbiosis as well as a satellite model for evolution in comparison to C. elegans. To explore gene expression differences across taxa, we sequenced the transcriptomes of single embryos of two Steinernema species and two Caenorhabditis species at 11 stages during embryonic development and found several interesting features. Our findings show that zygotic transcription initiates at different developmental stages in each species, with the Steinernema species initiating transcription earlier than Caenorhabditis. We found that ortholog expression conservation during development is higher at the later embryonic stages than at the earlier ones. The surprisingly higher conservation of orthologous gene expression in later embryonic stages strongly suggests a funnel-shaped model of embryonic developmental gene expression divergence in nematodes. This work provides novel insight into embryonic development across distantly related nematode species and demonstrates that the mechanisms controlling early development are more diverse than previously thought at the transcriptional level.

ALG-5 is a miRNA-associated Argonaute required for proper developmental timing in the Caenorhabditis elegans germline.

Caenorhabditis elegans contains 25 Argonautes, of which, ALG-1 and ALG-2 are known to primarily interact with miRNAs. ALG-5 belongs to the AGO subfamily of Argonautes that includes ALG-1 and ALG-2, but its role in small RNA pathways is unknown. We analyzed by high-throughput sequencing the small RNAs associated with ALG-5, ALG-1 and ALG-2, as well as changes in mRNA expression in alg-5, alg-1 and alg-2 mutants. We show that ALG-5 defines a distinct branch of the miRNA pathway affecting the expression of genes involved in immunity, defense, and development. In contrast to ALG-1 and ALG-2, which associate with most miRNAs and have general roles throughout development, ALG-5 interacts with only a small subset of miRNAs and is specifically expressed in the germline where it localizes alongside the piRNA and siRNA machinery at P granules. alg-5 is required for optimal fertility and mutations in alg-5 lead to a precocious transition from spermatogenesis to oogenesis. Our results provide a near-comprehensive analysis of miRNA-Argonaute interactions in C. elegans and reveal a new role for miRNAs in the germline.

Quantitative Analysis of Synthetic Cell Lineage Tracing Using Nuclease Barcoding.

Lineage tracing by the determination and mapping of progeny arising from single cells is an important approach enabling the elucidation of mechanisms underlying diverse biological processes ranging from development to disease. We developed a dynamic sequence-based barcode system for synthetic lineage tracing and have demonstrated its performance in C. elegans, a model organism whose lineage tree is well established. The strategy we use creates lineage trees based upon the introduction of synthetically controlled mutations into cells and the propagation of these mutations to daughter cells at each cell division. We analyzed this experimental proof of concept along with a corresponding simulation and analytical model to gain a deeper understanding of the coding capacity of the system. Our results provide specific bounds on the fidelity of lineage tracing using such approaches.

Linker histones: novel insights into structure-specific recognition of the nucleosome.

Linker histones (H1s) are a primary component of metazoan chromatin, fulfilling numerous functions, both in vitro and in vivo, including stabilizing the wrapping of DNA around the nucleosome, promoting folding and assembly of higher order chromatin structures, influencing nucleosome spacing on DNA, and regulating specific gene expression. However, many molecular details of how H1 binds to nucleosomes and recognizes unique structural features on the nucleosome surface remain undefined. Numerous, confounding studies are complicated not only by experimental limitations, but the use of different linker histone isoforms and nucleosome constructions. This review summarizes the decades of research that has resulted in several models of H1 association with nucleosomes, with a focus on recent advances that suggest multiple modes of H1 interaction in chromatin, while highlighting the remaining questions.

CeNDR, the Caenorhabditis elegans natural diversity resource.

Studies in model organisms have yielded considerable insights into the etiology of disease and our understanding of evolutionary processes. Caenorhabditis elegans is among the most powerful model organisms used to understand biology. However, C. elegans is not used as extensively as other model organisms to investigate how natural variation shapes traits, especially through the use of genome-wide association (GWA) analyses. Here, we introduce a new platform, the C. elegans Natural Diversity Resource (CeNDR) to enable statistical genetics and genomics studies of C. elegans and to connect the results to human disease. CeNDR provides the research community with wild strains, genome-wide sequence and variant data for every strain, and a GWA mapping portal for studying natural variation in C. elegans Additionally, researchers outside of the C. elegans community can benefit from public mappings and integrated tools for comparative analyses. CeNDR uses several databases that are continually updated through the addition of new strains, sequencing data, and association mapping results. The CeNDR data are accessible through a freely available web portal located at http://www.elegansvariation.org or through an application programming interface.